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A general model is presented for coupling of high-Q whispering-gallery modes in optical micro- 
sphere resonators with coupler devices possessing discrete and continuous spectrum of propagating 
00 ' modes. By contrast to conventional high-Q optical cavities, in microspheres independence of high 

On ' intrinsic quality-factor and controllable parameters of coupling via evanescent field offer variety of 

regimes earlier available in RF devices. The theory is applied to the earlier-reported data on differ- 
ent types of couplers to microsphere resonators and complemented by experimental demonstration 
of enhanced coupling efficiency (about 80%) and variable loading regimes with Q > 10* fused silica 
microspheres. 
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I. INTRODUCTION 



High-Q optical microsphere resonators currently attract growing interest in experimental cavity QED , measure- 
ment science, [||J^ frequency stabilization and other photonics applications Stemming from extensive studies of 
Mie resonances in microdroplets of aerosols |^ (observed via elastic and inelastic scattering of free-space beams), further 
studies of laboratory-fabricated solid-state microspheres focus on the properties and applications of highly-confined 
whispering-gallery (WG) modes. The modes of this type possess negligible electrodynamically-defined radiative losses 
(the corresponding radiative quality factors Qrad > 10^" and higher), and are not accessible by free-space beams 
(-H and therefore require employment of near-field coupler devices. By present moment, in addition to the well known 
O ^ prism coupler with frustrated total internal reflection 1^,1^, demonstrated coupler devices include sidepolished fiber 
'— ' \ coupler P jTl| , p^ and fiber taper . The principle of all these devices is based on providing efficient energy transfer 
^ • to the resonant circular total-internal-reflection-guided wave in the resonator (representing WG mode) through the 
] evanescent fleld of a guided wave or a TIR spot in the coupler. 

, It is evident a priori that efficient coupling maybe expected upon fulfillment of two main conditions: l)phase 
' synchronism and 2)significant overlap of the two waves modelling the WG mode and coupler mode respectively. 
; Although reasonable coupling efficiency has been demonstrated with three types of devices (up to few tens of percents 
of input beam energy absorbed in a mode upon resonance), no systematic theoretical approach has been developed 
^5 1 to quantify the performance of coupler devices. It still remained unclear whether it is possible at all and what are 
conditions to provide complete exchange of energy between a propagating mode in a coupler device and the given 
whispering gallery mode in high-Q microsphere. Answers to these questions are of critical importance for photonics 
^ . applications and also for the proposed cavity QED experiments with microspheres. 
. ^ ' In this paper, we present a general approach to describe the near-field coupling of high-Q whispering-gallery mode to 
a propagating mode in dielectric prism, slab or waveguide structure. Theoretical results present a complete description 
(-H ' and give a recipe to obtain optimal coupling with existing devices. We emphasize the importance of the introduced 
O ^ loading quality-factor parameter Qc and its relation with the intrinsic Qp-factor of WG modes as crucial condition to 
obtain optimal coupling. Theoretical consideration is complemented by experimental tests of variable loading regimes 
. ^ and demonstration of improved coupling efficiency with prism coupler. 

, II. GENERAL CONSIDERATIONS 



Let us examine excitation of a single high-Q whispering gallery mode with high quality-factor by (N) travelling 
modes in an evanescent wave coupler. This coupler can either have infinite number of spatial modes (N — oo, as in 
prism coupler PJlO[| and slab) or only one mode {N = 1, as in tapered fiber and integrated channel waveguide). 

We shall start with simple description of the system using lumped parameters and quasigeometrical approximation. 

Let Ao{t) be the amplitude of a circulating mode of total internal reflection in the resonator (see Fig.l) to model 
the whispering-gallery mode. Let the pump power be distributed in the coupler between its modes so that Bk{t) is 
the amplitude of mode k {1 < k < N) and B^"-{t) is the slow varying amplitude with ^ = |i?™(t)P equal to 

total pump power. Let us assume for simplicity that coupling between different modes is absent without the resonator. 
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Assuming that the couphng zone is much smaller than the diameter D of the resonator, we can introduce local 
real amplitude coefhcients of transmittance Tk to describe the coupling of the resonator with all modes of the coupler 
(either guided or leakage ones) and the internal reflectivity coefficient R. We shall denote arrays of transmittance 
coefficients and amplitudes as vectors T and B respectively. If the quality-factor of the resonator mode is high enough, 
then a single circle gives only small contribution to the mode buildup and therefore 1 — i? ^ 1. In this case (neglecting 
for simplicity absorption and scattering losses in the coupler Rf. = 1 — T|) we obtain 

R = l[Rk = l-Y.Ti/2 = l-Ty2 (1) 

The equation for the mode of the resonator will be the following: 

Ao{t) ^i'^TkB'jrit) + RAa{t - tq) exp [127^^1/ X - aL/2] 
Bl^\t)^RuBT{t)+iTMt). (2) 

where tq — UsL/c is the circulation time for the mode travelling inside the sphere, L ~ 27ra is approximately equal 
to the circumference of the sphere, A is the wavelength, ris is the refraction index, c is the speed of light and a is the 
linear attenuation in the resonator caused by scattering, absorption and radiation. 

In the above representation the microsphere is equivalent to a ring resonator formed by mirrors with transmittances 
Tk and filled with lossy medium, or, in case of single-mode coupling, as pointed in to a Fabry-Perot resonator of 
the length L/2 with totally reflecting rear mirror. 

If propagation losses are small, then near the resonance frequency ljq — 27rc/Ao, UgL = toAq, where m is integer, 
by expanding Ao(i — tq) = Ao(t) — TodAo/dt from (0) we obtain: 



dt 
where 



+iSc + S„ + iAu)Ao^iCB''^, (3) 
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We introduce here another important coefficient: 

r = (5) 

This coefficient F <1 describes mode matching and shows how close the field in the couplers matches the near field 
of resonator mode. 

The term 5q originates from intrinsic quality factor Qo = 27rns/aA while 5c describes loading i.e. mode energy 
escape to all modes of the coupler. Hereafter we shall mark all values associated with coupler by index 'c' and values 
associated with microsphere by index 's'. 

Equation (^) is a classical equation for the amplitude of the resonator pumped by harmonical field. 

As will be shown below, coefficients Tk can be calculated as normalized overlap integrals of the fields of the 
microsphere mode and modes of the coupler. The difference from Fabry-Perot resonator is that for the microsphere, 
coefficients Tk are not fixed parameters but instead, strongly depend on geometry of coupling (e.g. exponentially 
on the value of the gap between microsphere and coupler) and are therefore in hands of experimentalist. As we 
already emphasized in ||9|] , it is the controllable relation between 5q and 5c that defines coupling efflciency upon given 
configuration (accounting both for mode overlap and synchronism and optimized loading to provide energy exchange 
between resonator and coupler). Stationary solution for (^) has the typical form: 

^ _ i25cB''^ r _ iTB'"' [25^ 
5o + 5c + iAui T 5a + 5c + iAcj V tq 

Field amplitude in the resonator will be maximal at 5c = 5o (intrinsic quality- factor equals the loaded Q). The output 
stationary amplitudes are 

25 r T 

•QOut "Q^^^ B^^ 

6a + 5c + iAuj T 



and total output intensity in this case has lorentzian shape: 
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It can be easily seen from this equation that the output signal can be considered as the result of interference of the 
input and the "re-emission" from the resonator. Note that mode distribution of the second term in (resonator 
mode emission pattern) does not depend on the input distribution. 

The most important case of (0) is the regime of ideal matching (F = 1), obtained with B^Y^j™ = T/T when the 
fraction of the input power fed into the resonator mode is maximal. (Single- mode coupler is always "mode-matched" .) 
In this case, provided Sc = Sq, output intensity turns to zero i.e. the entire input power is lost inside the resonator. 
This regime is usually called critical coupling. Sometimes coupling is characterized by the fractional depth K of the 
resonance dip in intensity transmittance observed upon varying the frequency of exciting wave in the coupler; from 
K can be expressed as follows 

(Qo + Qc)2 Qo + Qc 
1 2^0 1 1 

n = — + — = 7r + 7r (9) 

In case of critical coupling K — 1 (100%). In case of nonideal matching, critical coupling may be observed until 
2r^ > 1 (partial matching) if the output is mode- filtered to pick up only part of the coupler modes. In this case 
leakage into other modes may be considered as additional internal losses, and critical coupling is obtained with lower 
loaded quality factor when Sc = 5o/{2r'^ — 1). If 5c ^ <5o (overcoupling) then for matched coupling the output wave in 
resonance has the sign opposite to that out of resonance i.e. the resonator shifts phase by tt. It is appropriate to note 
here that in traditional high-Q optical resonators comprised of mirrors, the quality-factor is limited by the mirrors' 
finesse i.e. by loading. With microspheres, the situation is opposite, and the primary role belongs to the intrinsic 
quality-factor. 



III. DIRECTIONAL COUPLER APPROACH 



The goal of this section is to determine parameters of the coupler-resonator system from the electrodynamical point 
of view. In the recent paper |l^] by D.R.Rowland and J.D.Love, famous for their popular book on the theory of 
optical waveguides ||l^ , the problem of coupling with whispering-gallery modes is addressed on the basis of the model 
of distributed coupling between a travelling surface mode in cylindrical resonator and a given mode in a planar (slab) 
waveguide. In this approach, the coupling problem leads to the necessity to solve a system of differential equations, 
which in our designations looks as follows: 

^ = zA/3o(2)Ao -I- iCk{z) exp - Po)z]Bk 



dz 
dBk 
dz 



iAPk{z)Bk + iCkiz) exp [-i{f3k - I3q)z]Ao, (10) 



Coefficients A(3o{z) and A/3fc(z) (describing perturbation of wave numbers /3o and Pk of modes of the resonator and 
the coupler) and distributed coupling coefficients Ck can be calculated explicitly as field cross-section integrals (see 
| p5[ and references therein). 

A/?o = J \eo\^ds; A/3.(.) = J [e.^ds; 

c s 

.2 



-^j{nl~ l)e*eods j [n] - l)e*ekds (11) 
c s 

Here Bq and are equivalent waveguide modes of the resonator and of the coupler respectively, normalized vs. power; 
the integration is done over cross-sections. Indexes S and C denote that the integration is done inside the microsphere 
and coupler respectively. In principle, conservation of energy requires that the two integrals in expression for be 
equal and this is frequently postulated. However, in common approximation that we also use here the above equality is 
secured only for phase-matched or identical waveguides, while in the opposite case the dependence of the two integrals 
on the gap is different. Nevertheless, to provide efficient coupling, this equality must be satisfied. 
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Parameters (11) are nonzero only in the coupling zone. It may seem that the coupler transmission matrix (CTM) 
and, subsequently, the above introduced lumped Tk coefficients can be found from equations (10). However, analytical 
derivation of the output field amplitudes cannot be found from (^0|) with exception of few simple cases. It was perhaps 
due to this fact that the authors of |jl^ presented only numerical solution for their particular case. Moreover, in general 
case CTM is a complex 2x2 matrix and cannot be characterized by one real parameter. 

Fortunately, the situation is more favorable for optical microsphere resonators with high loaded quality-factor 
Qc — uj/2Sc, when Tk ^ 1. Indeed, from it follows that — m/Q. In a fused silica resonator with the diameter 
140/i-TO (to ~ 1000) and heavily loaded Q ~ 10^ (intrinsic Q - factor can be of the order of 10^" in this case) T ~ 1%. 
In practice T is usually of the order of lO"*^. It means that the field amplitude Aq changes insignificantly over the 
coupling zone and can therefore be assumed constant in the second equation of (|o|), and the stationary amplitude 
Aq ^ Bf;. Therefore an approximate solution can be obtained: 



A°''' = exp 



^Bl 



iTkAo, 



iTkBl 



(12) 



where 



Tk= Ck exp [i{f3o - f3k)z]dz; 



(13) 



Equations (^2|) are practically identical to (||) if Aq is closed into a ring. In the second equation of ( |l^ ) we neglected 
small second-order terms while, however, keeping them in the first equation as they describe the coupler-induced shift 
in resonant frequency and the reduction of Q by loading. 



To 



$ _ m ini ^ 1) 
Sttto 



\eo\'^dv 



2to 1287r2To 



2^.eo exp [i(/3o - I3k)z\dv 



(14) 



IV. VARIATIONAL APPROACH 



Directional coupler approach can be easily generalized for multimode coupler. However expressions for the coupling 
parameters better suited for couplers with dense mode spectrum can be found in a more rigorous way directly from 
Maxwell equations using variational methods. Electric field in the resonator perturbed by coupler may be written in 
the form: 

E,(r,t) = exp [-iut] ^ ij(i)e,(r), (15) 

j 

where are orthonormalized eigenmodes of the unperturbed lossless resonator without coupler 

£s%ie*2dw = (^jij2 (16) 



i^ji,j2 is Kronecker symbol here). Rigorously speaking this normalization meets some difficulties for open dielectric 
resonators with finite radiative quality- factor Qrad [0- In our consideration, however, we can avoid them by assuming 
eigenfrequences of interest cuj to be purely real. For this we neglect the imaginary part that describes radiation losses 
and choose as the integration volume the sphere with a diameter much less than QradA/7r. Amplitudes Aj are slowly 
varying and differ from circulating amplitude Aj introduced before only in terms of normalization. One can easily 
see that 

pl = ^f [e,h*]ds - -. (17) 
|ej|2 AttJ ' ^ ^' Tj 
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The equation for the field in the coupled sphere will have the form: 

VxVxE,+ (^-^ + ^^ + z— -^j^- (18) 

where the second term in brackets is additional polarization due to presence of the coupler, the third one describes 
damping associated with intrinsic losses in the resonator, and the right part is the polarization caused by the pump 
wave. Dielectric susceptibilities es|c(r) are equal to inside and unity outside the spherical resonator and the 
coupler correspondingly. Substituting ( p^ into ( p^ ) and multiplying this equation by Gq, after integration over the 
entire volume and omitting small terms we obtain: 

^ + Aoido + iAiu') = exp [iut] I E^Bldv, (19) 



where Aoj' = lu'q — lu and 



in 



co'o=ujo-^ I {nl-l)\eo\^dv. (20) 



is the new resonance frequency shifted due to the coupler, in total agreement with (^^. Let us express the field in 
the coupler as expansion in travelling modes in z direction: 

Ec(r, t) = exp [~iLut] J Bp{z, t)e^(r) exp [il3z]d(3, (21) 

Guided localized modes of the coupler in this description can also be easily taken into account if we choose Bp as 

= ^Sfc<5(/3-/?fc)+B^ (22) 

k 

The coupler modes are normalized in such a way that 

^ J[epi,h}j^],ds = S{Pi-P2) (23) 

(here S{Pi — P2) is delta-function and h is the magnetic field corresponding to the mode). Integration is performed 
over the cross-section orthogonal to z-axis. Amplitudes B/3 (slowly varying with z and t) describe distribution of the 
pump wave in coupler modes. Substituting ( pT| ) into the wave equation: 

we obtain 

- ^^%^^/^) -P i^^-W - ^^^^^^E. exp [^u.t] (25) 

The second term in brackets (|2^ ) determines the change of the wavenumber (phase velocity) for the given mode in 
the coupling zone. Taking vector product of this equation with and integrating over the cross-section, we obtain 
formal solutions for slowly varying amplitudes: 

Bp = Bf exp [iA/3z] + j exp [iut - i(3z' + iA/3(z - z')] j {nl - l)[¥.'^,h'^]M dz' 

-00 c 

A/3(z) = ^ j [e^, h^l.ds, (26) 

s 

Substituting (|l|) into ( p^ ) using ( p6| ) and omitting A/3, we finally obtain the following equation for the amplitude 
of the mode in the resonator: 
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d-A-o _ ^ . A /\ J iuj(n'^ — 1) 
+ (^o + 5c + iAlu')Ao = '- 

at OTT 



B™ / eoe*pexp{—iPz)dvd(3, 



(27) 



and 



647r2c 



S -oo C 



1287r2 



ege^ exp [—i^zjdv 



dp 



(28) 



in natural agreement with ([T^). Total agreement with (^|-0) becomes apparent if we put 



^{nj - 1) 
Btt 



eoe^ exp [z(/3o - P)z]dv 



(29) 



For high-Q WG modes f3o ~ m/a, and as the field drops outside the resonator approximately as exp [— 7r] (7^ ~ 
fc^(n^ — 1)), the dependence of eo on z can be approximated as follows: e ~ e(z = 0) exp [—72^/20]. If the coupler is 
straight in z direction (as in most demonstrated couplers to date), we obtain: 



Tb 



Lj{ni - 1) l2na 
8^ 



exp [—(/3a — m)'^ /2^a]eoe%ds 



(30) 



c 



V. APPLICATION TO DEMONSTRATED COUPLERS 



Let us now use the developed approach for the analysis of coupling of whispering-gallery modes with optical fiber. 
As soon as according to (refkcoupl), possibility of efficient coupling critically depends on the value of the loading 
quality- factor Qc and its relation with the intrinsic Q, in this section we shall focus on calculation of Qc and discuss 
briefly methods to achieve phase synchronism and mode matching with different couplers. 

To date, two types of optical fiber coupler to WG modes in microsphere were demonstrated. The first one is the 
eroded fiber coupler where evanescent field of a propagating waveguide mode becomes accessible due to 

partial removal of the cladding in a bent section of the fiber. The recently demonstrated second type of the fiber 
coupler is based on the stretched section of a single mode fiber employing the mode conversion of the initial guided 
wave into waveguide modes of cladding tapered to the diameter of few microns . 

The most interesting type of strong confinement modes of the sphere - TEuq with radius a (where radial index q is 
small) and the HEn mode in the fiber of radius b can be approximated as follows |r^ , [l^ : 



'^V / Jo(w)/Jo(ryfe) ifp<6 



where 



-^b^/nc \ exp [-^{p - 6)] if p > 6 



r?6 ~ 2.405 exp -— ^ — ~ 2.405 1 - ' 



(32) 



267 y ■ V 267 



7 ~ V(^+l/2)Va'-fc' - fcv/«2 - 1 (33) 
Using ( |l^ and several approximations, we can now calculate 

Qc^ 9(n2- l) ''''P [2kd^n^ - I + {t - Paf/ia], (34) 
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where d is the gap between the resonator and the fiber and Us = ric = n. To obtain optimal coupHng, one has to 
require matching the propagation constants in the argument of the second exponent of (|3^ ) {£ = 13a) as in [p^ . In 
this case, using approximations for eigenfrequences in the resonator, one can obtain optimal radius of the fiber and 
the loaded Q. 

2.3a „ / aA2 ^ 



0.51 



Q, ~ 102 (^) exp [2Vn^ - Ikd] (35) 

VA/ 4g — 1 

Using the above expressions, we can try to compare our calculations with the experimental data reported in [T^ ] 
for A — 1.55 /im, a = 85 /im and b = 1.7 ^m. The measured Q was 2 x 10^ with K = 72%. Using we can obtain 
Qo — 8.5 X 10^ and Qc = 2.6 x 10^. Calculations with ( ^ ) give Qc = 2.5 x 10^ - in agreement with the experiment. 

It is appropriate to note here that in principle, as follows from (|3^), the minimum of Qc does not correspond to 
phasematching (/3 = £/a) and is shifted to smaller b. This minimum is also not very sharp (~ 2/^/1 - several percents 
of However this case deserves special consideration, because for smaller b the approximations we use here will 
give larger (more than 10%) error. It is also important that the loaded Q increases very quickly with the size of 
the resonator (as £^/^), and in this way the range of possible applications of such coupler becomes restricted. Even 
for very small fuscd-silica spheres, optimal radius of most common silica fiber does not correspond to single-mode 
operation implying further technical complications in using this type of coupler. 

The conclusions of our theory also correlate with the data on limited efficiency of sidepolished optical fiber couplers 
P, p]Jl^ . Indeed, with the typical monomode fibers having core index equal or smaller than that of the spheres (made 
from polystyrene or silica correspondingly in 11 1^, with small microspheres one cannot satisfy phasematching 



because of relatively large diameter of standard cores, and with larger spheres coupling coefficient is too small (coupling 
Q-factor too high) to provide efficient power insertion into the resonator. 

Efficient coupling (tens of %) with high-Q microspheres has been demonstrated with the planar (slab) waveguides 
p9[ . This type of the coupler provides additional freedom compared to fiber waveguides because it allows free 
manipulation of the two-dimensional optical beams. Optimal width for ££q mode is g = 2a/V2l. In the meantime 
requirement of phase matching for efficient coupling implies optimization of the slab waveguide thickness /. Using 
the same approach as above, we obtain: 



^ ^, 8A3r.,(7;i|~T+n,)(n^-l)3/^ a/3 

Qc ^ _ ly JT [2fcrfV"i + l3a) /ja] (36) 

It is appropriate to note here that either fibers or planar waveguides can also effectively excite modes with £ ^ to if 
the wavevector is inclined to the "equator" plane of the sphere (symmetry plane of the residual ellipticity) by angle 
arccos(TO/£). This conclusion becomes evident if we remember that the mode with Imq is equivalent to a precessing 
inclined fundamental ££q mode . 

Prism coupler has been analyzed in our previous papers ]9|JTo|] together with precession approach to the description of 
WG modes and theoretical and experimental investigation of the far field patterns. By contrast to waveguide couplers, 
where practical realization of high efficiency implies either precise engineering of the waveguide parameters, or the step- 
by-step search of the optimal contact point to the fiber taper, prism coupler allows systematic procedure of coupling 
optimization by manipulating the external beams. The two steps to achieve efficient coupling are l)adjustment of the 
incidence angle $ of the input Gaussian beam inside the coupler prism and 2) adjustment of the angular size of the 
beam A$ and A8 to provide mode matching with far field of the WG mode in the prism {F factor). 



£ . n x/nl — 1 . T ns + \/n'i — 1 



The loading Q with the prism coupler is as follows: 

Q-- ^ /.Sjl \ d -p[2fcdV^] (38) 



Fig. 2 summarizes the calculations of the loading quality-factor for different types of couplers (with optimized 
parameters) in form of the plots of Qc under zero gap d = 0. The results in Fig. 2 allow to quickly evaluate possibility 
to achieve critical coupling with the given size and intrinsic Q of the sphere, along the lines summarized in Sec. 2 (||). 
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In our experiments employing high-Q WG microsphere resonators, we used prism coupler in most cases and believe 
that it remains the most flexible device as it provides ability of fine adjustment of both phase synchronism and mode 
matching via convenient manipulation of the apertures and incidence angles of free beams. Also, as seen in Fig. 2, it 
provides a significant margin to obtain critical coupling with the spheres of various size and intrinsic Q. As a result, we 
routinely obtained coupling efhciencies to silica microspheres about 30% with standard right-angle BK7 glass prisms 
limited by restricted mode overlap due to input refraction distortions of the symmetrical Gaussian beams. Use of 
cylindrical optics or higher refraction prisms to eliminate the mode mismatch (F —^ 1) can significantly improve 
coupling efhciency and approach full exchange of energy under critical regime. (About 80% coupling efficiency is 
demonstrated further in the experimental Sec. 7). 

To conclude this section, let us also note here that the critical coupling (that is characterized by maximal absorption 
of the input power in the resonator) is in fact useless for such applications as cavity QED or quantum-nondemolition 
experiments, because no recirculated power escapes the cavity mode. To be able to utilize the recirculated light, one 
has to provide the inequality Qc ^ Qq (strong overcoupling) . In other words, the intrinsic quality- fact or has to be 
high enough to provide reserve for sufficient loading by the optimal coupler. 



VI. PRISM COUPLER: EXPERIMENTAL EFFICIENCY AND VARIABLE LOADING REGIMES 

In order to illustrate the results of our analysis, we performed measurements to characterize coupling efficiency 
of the prism coupler with high-Q fused silica microspheres. As in our previous experiments, we used microspheres 
fabricated of high-purity silica preforms by fusion in a hydrogen-oxygen microburner (see description of the technique 
in pc|]). The focus of present experiment was to obtain enhanced coupling by maximizing the mode matching 



along with the lines briefly described in Sec. 2. In our experiment, to diminish astigmatic distortions of the input 
gaussian beam at the entrance face of the coupler prism, we used equilateral prism of flint glass (SF-18, refraction 
index n= 1.72). As usual, the input beam (a symmetrical gaussian beam from a single-mode piezo-tunable He-Ne 
laser) was focused onto the inner surface of the prism, at the proximity point with the microsphere. The angle of 
incidence and the cross-section of the input beam were then optimized to obtain maximum response of a chosen WG 
mode. Initial alignment was done on the basis of direct observation of resonance interference in the far field, with 
the frequency of the laser slowly swept across the resonance frequency of the mode. With the given choice of prism 
material, optimal angle of incidence for excitation of whispering-gallery modes g ~ 1 (close to critical angle of total 
internal refraction at the silica-glass interface) was approximately equal to 60 degrees so that astigmatic distortions 
of the input beam at the entrance face of the prism were minimized. 

After preliminary alignment, the coupling efficiency was further maximized on the basis of direct observation of the 
resonance absorption dip: full intensity of the beam after the coupler prism was registered by linear photodetector and 
monitored on digital oscilloscope. Results obtained with a TMuq mode (possessing strongest confinement of the field 
in meridional direction) are presented in Fig. 3 in form of the resonance curves observed upon successively decreasing 
coupling (stepwise increasing gap). Fig. 3 illustrates good agreement of theory with experiment: indeed, resonance 
transmission decreases with loading until the quality- factor becomes twice smaller than the intrinsic Qo! after that, 
intensity contrast of the resonance decreases. Fig. 4 presents explicitly the plot of the fitted experimental intensity dip 
versus the loaded quality- factor Q = QcQo/{Qc + Qo), which yields satisfactory agreement with parabofic prediction 
from the generalized expression (0). Maximal contrast of the resonance obtained in our experiment was ~ 0.79 
(the "deepest" curve in Fig. 3). 



VII. CONCLUSION 



We have presented a general approach to describe the near-field coupling of high-Q whispering-gallery modes in 
optical microsphere resonators to guided or free-space travelling waves in coupler devices with continuous and discrete 
spectrum of propagating modes. 

A convenient formalism of the loaded quality-factor to describe the energy exchange between coupler modes and 
the resonator provides a quick algorithm to determine the efficiency of the given type of the coupler, under given 
value of the intrinsic quality-factor of WG modes. 

Variable relation between the intrinsic Q-factor and loading losses (described by Qc) through energy escape to 
coupler modes is a distinctive new property of whispering-gallery resonators compared to conventional Fabry-Perot 
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cavities: the latter are characterized by fixed couphng through the reflectivities of comprising low-loss mirrors. This 
unique ability to control the Q and coupling via WG mode evanescent field allows to obtain new regimes in the 
devices, analogous to those available in lumped-element RF and microwave engineering. 

Theoretical estimates on the basis of the suggested theory are in good agreement with the reported data on the 
efficiency of different coupler devices including tapered, sidepolished fiber and slab waveguide. 

Original experimental results include direct demonstration of variable loading and enhanced efficiency (up to about 
80% ) in prism coupler. Ease of control of phase synchronism and mode overlap between coupler and microsphere 
mode by adjusting the input beam parameters make the prism coupler versatile and efficient for various applications 
of high-Q microsphere resonators. 

In conclusion, let us note that the near-field coupling may be not a unique method to efficiently excite highly confined 
whispering-gallery modes in microspheres. Simple estimates show that for example, recent advances in optical fiber 
grating fabrication methods ||2l| may allow to " imprint" a Bragg-type critical coupler for high-Q WG modes directly 
on a sphere made of low-loss germanosilicate glass. This configuration might be of special interest for atomic cavity- 
QED experiments, where presence of bulky external couplers may destroy the field symmetry, complicate laser cooling 
of atoms etc. 
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FIG. 1. Schematic of the excitation of whispering gallery modes in high-Q microsphere 
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Microsphere radius, alX 



FIG. 2. Efficiency of different couplers in terms of tlie loading quality-factor Qc witli optimized parameters at d = as 
function of spliere radius a; numerical results are obtained for TMui mode. Critical coupling is possible when the intrinsic 
quality-factor of WG mode Qo is larger than Qc (see (P)) 
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FIG. 3. Output intensity of the prism coupler observed under variable loading (successively increasing microsphere-prism 
gap). Fused silica sphere with diameter 270^m. 
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